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We report on the detection with the Infrared Space Observatory (ISO) of a 
molecular band at 57.5 /im (174 cm -1 ) in carbon-rich evolved stars and in Sgr B2. 
Taking into account the chemistry of these objects the most likelihood carrier is 
a carbon chain. We tentatively assign the band to the 1/5 bending mode of C4 for 
which a wavenumber of 170-172.4 cm -1 has been derived in matrix experiments 
(Withey et al. 1991). An alternate carrier might be C4H, although the frequency 
of its lowest energy vibrational bending mode, is 7 , is poorly known (130-226 
cm -1 ). If the carrier is C4, the derived maximum abundance is nearly similar to 
that found for C3 in the interstellar and circumstellar media by Cernicharo et al. 
(2000). Hence, tetra-atomic carbon could be one of the most abundant carbon 
chain molecules in these media. 

Subject headings: infrared: general — ISM: individual (Sgr B2) — ISM: lines 
and bands — ISM: molecules — line: identification — stars: individual (IRC+10216, 
CRL2688, CRL618, NGC7027) 
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1. Introduction 

Long carbon chain radicals have been detected in interstellar and circumstellar clouds 
through their pure rotational spectrum at radio wavelengths (C5H, CeH, C7H, CsH, H2C3, 
H 2 C 4 : Cernicharo et al. 1986, 1987, Guelin et al. 1987, 1997, Cernicharo and Guelin 1996). 
Polyynes such as C 2 H 2 , C 4 H 2 , C 6 H 2 and single aromatic species such as benzene, have been 
detected in the protoplanetary nebula CRL618 (Cernicharo et al. 2001a&b) supporting the 
idea that much complex carbon-rich molecules could be formed in space. However, the 
mechanisms allowing the growth of carbon-rich molecules are still poorly known and the full 
set of molecules that could participate in the chemical reactions leading to the formation of 
large complex carbon-rich species has yet to be identified. Among these "building blocks" 
the C n chains have attracted the interest of laboratory spectroscopists and astronomers. 
C3 has been observed at optical wavelengths in the atmospheres of cool stars (see, e.g., 
Zuckerman et al. 1976), and it has been identified in the envelope of IRC+10216 through its 
z/ 3 antisymmetric stretching mode in the mid-infrared by Hinkle et al. (1988) who derived 
X(C 3 )~1CT 6 . This key species has recently been found in the ISM (Cernicharo et al. 1996, 
1997a, 2000, Giesen et al. 2001) through its z/ 2 low frequency bending mode (Schmuttenamer 
et al. 1990) and in the diffuse interstellar clouds through its electronic transition near 4052 
A (Maier et al. 2001; Roueff et al. 2002). It has been suggested that C3 could be involved 
in the formation of the diffuse interstellar bands (Douglas, 1977; Clegg and Lambert 1982). 
ISM observations in Cernicharo et al. (2000) put a lower limit of 10 -8 to the C3 abundance. 
In addition, C5 has been detected in the C-rich star IRC+10216 with an abundance ratio 
C3/C5 ~10 (Bernath et al. 1989). Less is known about numbered carbon chains, but C 4 
could well be a very abundant molecule both in C-rich evolved stars and in the ISM. 

In this Letter we report the detection of a blend of lines at 57.5 /xm that we tentatively 
assign to the z/ 5 bending mode of C 4 . The feature has been observed in Sgr B2, IRC+10216, 
CRL618, CRL2688, and NGC7027. We discuss the possibility that the carrier could be 
another different carbon chain : C 4 H, although photodetachement experiments suggest that 
its lowest energy bending mode could be at higher frequencies (Taylor et al. 1998). 

2. Observations and Results 

We have searched for the emission/absorption of the low energy bending modes of 
polyatomic molecules toward interstellar and circumstellar clouds using the Long Wavelength 
spectrometer (LWS; Clegg et al. 1996, Swinyard et al. 1996) on board ISO 1 (Kessler et al. 
1996). Our main target for this search has been the C-rich evolved star IRC+10216. 
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The LWS grating data of IRC+10216 were analyzed and modeled by Cernicharo et al. 
(1996). The data have been reprocessed following the last pipeline product (TDT 19800158). 
No significant differences have been found with the already published analysis. After identi- 
fying all the HCN (including several vibrational states) and CO lines, a conspicuous feature 
composed by a blend of several lines (see Figure 1, top panel) was identified between 56 
and 58 /zm. These lines could not be assigned to any known species and no other similar 
patterns were found in the ISO far-infrared spectrum of IRC+10216. We have checked for 
all pure rotational lines of light species that could be potentially abundant without success. 
Hence, the only possible explanation is that the feature corresponds to a bending mode of 
an abundant heavy molecule, and most likely to a carbon chain. We searched for the same 
feature toward other objects selected from the same criteria of our previous searches of C3, 
C 4 H 2 , C 6 H 2 , and benzene (Cernicharo et al. 1996, 1997a, 2000, 2001a,b), i.e., IRC+10216, 
CRL618, CRL2688, NGC7027, and Sgr B2 (among other sources). 

For Sgr B2 M we have used the LWS grating spectra obtained during revolutions 287 and 
494 (TDTs 28701401, 287012130, 287012131, and 49400302). Some of these data have been 
already presented (Cernicharo et al. 1997b, Goicoechea & Cernicharo 2001). All together 
represent a total of 32 individual scans. The feature at 57.5 /zm is clearly present. There are 
many other LWS grating observations of Sgr B2 with a reduced number of scans. We have 
inspected all of them and, again, the feature appears in each individual scan. Sgr B2 N was 
observed in the course of a fast mapping of the Sgr B2 star forming region with the LWS 
spectrometer and only four scans are available (see Cernicharo et al. 1997b). 

For the Planetary Nebula NGC7027 we have used all the LWS grating data taken by ISO 
(see Herpin et al. 2002 for details). The data for the Proto-Planetary Nebula (PPN) CRL618 
(TDTs 68800302 and 68800450) are those published by Herpin and Cernicharo (2000). For 
CRL2688, a young PPN, the TDT 02101504 data have been used. A previous analysis of 
the ISO/LWS spectrum of this source has been presented by Cox et al. (1996). 

All the spectra have been analyzed using the ISO spectrometers data reduction package 
ISAP 2 . The data products correspond to version 10 of the pipeline. A polynomial baseline 
had to be removed from the final spectrum in each source. For some objects this baseline 
removal is critical to correctly extract the shape and intensity of the absorption features. 
Due to the large number of lines in Sgr B2 M and N spectra, the determination of the 
continuum level is not obvious. For sources with emission lines, although baseline removal 
is less critical, the global behavior of detector SW2 does not follow the expected shape for 
continuum emission and a polynomial of fourth degree has been removed after blanking the 
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wavelength intervals where features above 5 sigma were present. 

Figure 1 shows the 53-62 /im spectra of the selected sources. We have also analyzed the 
available observations for several other sources : O-rich stars, Orion, Sgr A, bright molecular 
sources, and even Mars to check for possible instrumental effects. The new infrared feature 
was not detected in any of these sources. We have also examined the spectral response of 
the LWS-SW2 filter in order to check for possible instrumental contributions to the blend 
of lines found around 57 /im but there is nothing at this wavelength (Garcia-Lario, private 
communication). Hence, we conclude that the blend of lines is real. The pure rotational 
lines of CO and HCN are indicated in the top panel of Figure 1. The model of Cernicharo 
et al. (1996) predicted a low flux for these high-J lines. Only in IRC+10216 some of the CO 
lines (J=43-42, 44-43 and 45-44) are marginally detected at a 3a level. 

The observed band consists of a blend of at least 4 spectral lines or sub-bands (at 57.0, 
57.6, 57.9 and 58.2 /tin). There is an additional feature separated from the band at 56.15 
/im that is detected in C-rich evolved stars but is less evident in Sgr B2 N and M. The 
J=47-46 line of CO could only make a modest contribution to this feature in C-rich evolved 
stars since other high-J CO lines are weak - even in IRC+10216. In Sgr B2 M this feature 
is blended with the absorption produced by NH 3 (J=9-8). The 56.15 /im feature could be 
one of the fine structure lines of [SI] at 56.31 /im ( 3 Po- 3 Pi) but its companion [SI] ( 3 Pi- 3 P2) 
line at 25.249 /tm is missing in all objects. The feature could also be the R(l) line of HD. 
However, it is unlikely to find such a high intensity for this line in evolved stars. We think, 
that if instrumental, the feature could arise from Uranus, the planet used as calibrator for 
the LWS/grating spectrometer. However, its intensity changes by more than a factor 10 
between NGC7027 and IRC+10216 or CRL618 (the continuum flux at 58 /tin is nearly the 
same for NGC7027 and CRL618, 1.15 and 0.7 10~ 16 W cm" 2 /!™- 1 respectively, and 4 and 
3 10~ 16 W cm- 2 fim- 1 for IRC+10216 and CRL2688 respectively). Therefore, the 56.2 /im 
feature in IRC+10216, CRL618, NGC7027 and CRL2688 cannot be explained alone by a 
possible instrumental contribution from the calibrator. Figure 1 shows that the 56.2 /tm line 
follows the intensity of the 57-58 /im features. In the following we will consider that all them 
belongs to the same molecular band. 

Two additional features at 53.7 and 54.4 /tm in the spectrum of IRC+10216 (labelled 



1 Based on observations with ISO, an ESA project with instruments funded by ESA Member States (espe- 
cially the PI countries: France, Germany, the Netherlands and the United Kingdom) and with participation 
of ISAS and NASA. 

2 The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instruments 
Teams and Data Centers. Contributing institutes are CESR, IAS, IPAC,MPE, RAL, and SRON. 
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Ul and U2 in the top panel of Figure 1), and another one at 52.89 jum (not shown), are 
also detected in NGC7027 but not in the other sources. These three lines remain, so far, 
unidentified. The contribution of OH (Goicoechea & Cernicharo 2002) and the maximum 
possible contribution from H 2 (Cernicharo et al. 2002, in preparation) to the spectra of 
Sgr B2 N and M are indicated in the panel corresponding to Sgr B2 M. The contribution of 
H2O to the spectrum of CRL618 at 58/im is negligible (Herpin and Cernicharo 2000). 

3. Discussion 

The fact that the 58 /zm feature is detected in all the selected C-rich post-AGB objects 
clearly points toward a carbon-rich molecule as the carrier of the band. Several molecules 
with four atoms have low frequency bending modes around 160-220 cm -1 . HC3N has a low 
energy bending mode, v 7) at 222.413 cm -1 (~45 /zm). Dyacetylene, C4H2, recently detected 
in CRL618 with a very large abundance (see Cernicharo et al. 2001a&b), has a low energy 
bending mode, u 9 , also at ~45 /iin. Another example of molecule with bending modes in 
this wavelength range is C 4 H. Pure rotational lines from the v 7 — 1,2 bending states of C 4 H 
have been detected in IRC+10216 by Guelin et al. (1987) and Yamamoto et al. (1987). 
The frequency of this mode is expected from ab initio calculations (Graf et al. 2001) to be 
around 178 cm -1 . Yamamoto et al. (1987) have estimated a frequency of 131 cm -1 from the 
1-doubling constant, q7, of the u 7 =l,2 vibrational levels. A 20% of error could be expected 
on this value (Mikami et al. 1989). Photoelectron spectra of C4H indicated, however, a 
larger frequency for this mode (226 cm -1 ; Taylor et al. 1998). The low lying 2 IT electronic 
state of C 4 H is at ~468 cur 1 above the ground (Taylor et al. 1998). Both transitions, 
the v 7 = 1 — and the electronic A 2 Tl — A 2 E are in the wavelength coverage of the Short 
Wavelength Spectrometer (SWS) on-board ISO. However, IRC+10216 is so bright in the 
mid-infrared that the expected range for these transitions is covered by strong ripples and 
not obvious features have been found so far (Cernicharo et al. 1999). The shape of the band 
in Figure 1 allows to rule out a x n — 1 E transition, i.e., the bending modes of HC3N, HC5N, 
C n H 2 with n even, C n with n odd, etc. These molecules have a X S ground electronic state 
and, at the limited resolution of the LWS grating spectrometer, a strong Q-branch with weak 
P and R-branches are expected for their bending modes. The HC 3 N v 7 and C 4 H 2 u 9 modes 
are present in the LWS spectrum of CRL618 (to be published elsewhere), and correspond to 
these expectations. 

Taking into account the large abundance of C3 found in Sgr B2 and IRC+10216, C 4 
could be a good candidate to explain the IR feature under investigation. Linear C 4 was 
first observed in the laboratory using electron spin resonance spectroscopy of carbon clusters 
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in matrices by Graham et al. (1976). They were able to derive the spin-coupling and 
spin-doubling constants, A and 7, for the ground vibrational state. C 4 was first detected 
and characterized in the gas phase by Heath and Saykally (1991) by observing several ro- 
vibrational lines of the is 3 antisymmetric stretching mode. Moazzen-Ahmadi et al. (1994) 
observed the same band and the Z/3+Z/5— z/5 hot band and determined the rotational constants, 
B and D, for both the ground and the z/ 5 states. The 1-doubling constant, q5, for the latter 
vibrational state was also derived leading these authors to estimate a frequency of 160±4 
cm -1 for the z/ 5 mode. The z/ 5 bending mode of C4 has been measured in matrix experiments 
at a wavelength of 58 jum (Withey et al. 1991), very close to that of the newly found IR 
feature. They observed two different lines separated by 2 cm -1 (170 and 172.4 cm -1 ) that 
they interpreted as two different trapping sites for C4 in the Ar matrix. Both values could 
correspond to the frequency of the z/5 mode of C4 but the authors adopted the value of 
172.4 cm -1 . This value agrees well with the series of peaks observed in the photoelectron 
spectroscopy of Arnold et al. (1991). The difference between the frequency of the new IR 
band and that derived from the matrix experiments is 2.6 cm -1 , which could be consistent 
with a frequency shift introduced by the Ar matrix of 1.5%. Ab initio calculations predict 
wavenumbers in the 175-214 cm -1 range (47-57 /im) and an infrared intensity of 46.2 km 
mol" 1 [A(z/ 5 = 1 - 0) ~ 0.25 s -1 ; Martin et al. 1991; Watts et al. 1992; Taylor and Martin, 
1996; see also Van Orden and Saykally, 1998]. Due to the complexity of the C4 electronic 
energy levels and to the limited accuracy level of ab initio calculations, this intensity value 
has to be taken into account as a crude estimation. 

The fact that the observed band (see Figure 1) contains several spectral features sug- 
gests a transition between terms of higher spin multiplicity, in particular toward a vibronic 
transition 2 U — 2 X or 3 n — 3 X, i.e., just the type of bands we can expect for C4H and 
C4 which have 2 £ and 3 S~ ground electronic states respectively. Consequently the lowest 
bending modes of these species, z/ 7 for C4H and z/5 for C4, will have 2 I1 and 3 n vibronic 
character respectively. An argument against C4H is the following, while pure rotational lines 
of its bending mode have been detected towards IRC+10216 (Guelin et al. 1987; Yamamoto 
et al. 1987) these lines are missing in all other objects reported here (except in CRL618 
where several weak lines from these states have been detected; Cernicharo et al. 2002, in 
preparation). Moreover, the abundance of C4H in Sgr B2 is rather low (Guelin, Mezaoui & 
Friberg, 1982), while C3 has a large abundance (Cernicharo, Goicoechea and Caux, 2000). 
Hence, it is tempting to assume that the observed features in post-AGB stars and in the 
ISM correspond to the z/5 mode of C4, although a definitive assignment will require better 
spectral resolution data and/or additional laboratory information on this species. The main 
problem in assigning the band to the z/ 5 mode of C 4 is the large spin-orbit interaction con- 
stant, A50, that should have the 3 n vibronic state in order to reproduce the observations. 



-7- 



As the ground electronic state of C 4 is 3 £, a very small value for A S o should be expected. 
However, if the molecule has a low lying excited electronic state with A#0, then the ground 
and the excited states will be strongly coupled through vibronic mixing and a large spin- 
orbit constant could be expected. This is the case for CCH which has a 2 E ground state 
and an A 2 I1 state 3600 cm -1 above. The 1^2 bending mode of the ground electronic state has 
Aso —0.3 cm -1 . For C4H, A$o is very large, ~3 cm -1 (Yamamoto et al. 1987), because its 
A 2 n state is only ~468 cur 1 above the 2 E ground state (see above). For C 4 the situation 
could be very similar to that of C 4 H as its first excited electronic state is 1 A and is only 2680 
cm -1 above the ground electronic 3 S state (Arnold et al. 1991; Xu et al. 1997). Last panel 
in figure 1 shows the computed band shape for As =4 cm" 1 at different temperatures. The 
similarity between the observed and computed shape is appealing. Nevertheless, we note 
that the adopted spin-orbit constant is particulary large. We have run several models with 
different Aso values and conclude that the shape of the band is only crudely reproduced for 
Aso >2 cm -1 . Another possibility that may explain the complexity of the band shape and 
allows smaller Aso values is that 13 CCCC and C 13 CCC isotopomers also make a remarkable 
contribution to the observed emission. This could also apply to C 4 H if finally this molecule 
is confirmed as the carrier. We note, however, that in this case we should also expect a very 
large column density of C4H, which looks incompatible with the upper limits found for this 
molecule in Sgr B2 (Guelin, Friberg and Mezaoui 1982). If the carrier is C4 then the column 
density will be particularly large in all observed sources. Assuming a vibrational excitation 
temperature of 50-100 K in AGB stars and PPNs, N(C 4 ) has to be close to 10 15 cm -2 , i.e., a 
typical abundance of a few 10~ 7 - 10 -6 . In Sgr B2 the computed column density (assuming 
a lower vibrational excitation temperature) is 1.2 10 15 cm -2 , i.e., similar to that derived by 
Cernicharo, Goicoechea and Caux (2000) for C 3 . The accuracy of the band intensity from 
ab initio calculations dominates the error in the estimated column densities. 

Future space-born platforms, such as the Herschel Space Observatory, will be equipped 
with heterodyne instruments with much better spectral resolution which should permit the 
detection of longer C n chains through their low lying bending modes in the far-IR. The 
detection of these molecules and the determination of their abundances will allow to stablish 
their role in the formation of the carriers of the unidentified infrared bands. Waiting for 
Herschel, important work has to be done in spectroscopy laboratories to fully characterize 
these molecular species. 

We thank Spanish DGES and PNIE for funding support under grant ESP98-1351-E 
and PANAYA2000-1784. JRG acknowledges Universidad Autonoma de Madrid for a pre- 
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Fig. 1. — Observed ISO/LWS spectra toward the sources discussed in the text. The wave- 
lengths of the CO lines J=43-42 up to J=49-48 and those of HCN J=56-55 up to 64-63 are 
indicated by arrows in the top panel. For IRC+10216, CRL2688, CRL618 and NGC7027 
the continuum level has been removed. For Sgr B2 N and Sgr B2 M the line over continuum 
flux ratio is shown. In all cases the continuum has been derived from a polynomial fit to the 
spectra after removing all emission/absorption features above 3 a. The bottom panel shows 
the expected band shape for a transition 3 n — 3 X with Asq=4 cm -1 , intrinsic linewidth of 
10 kms -1 and spectral resolution, A/AA, of 200. 



- 11 - 




